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We consider a porous medium equation with nonlocal diffusion effects given 
by an inverse fractional Laplacian operator. In a previous paper we have found 
mass-preserving, nonnegative weak solutions of the equation satisfying energy 
estimates. The equation is posed in the whole space IR n . Here we establish the 
large-time behaviour. We first find selfsimilar nonnegative solutions by solving 
an elliptic obstacle problem with fractional Laplacian for the pair pressure- 
density, which we call obstacle Barenblatt solutions. The theory for elliptic 
fractional problems with obstacles has been recently established. We then use 
■ entropy methods to show that the asymptotic behavior of general finite-mass 

solutions is described after renormalization by these special solutions, which 
represent a surprising variation of the Barenblatt profiles of the standard porous 
medium model. 
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1 Introduction 



In the paper [TT] we have introduced a model of nonlinear diffusion with nonlocal 
effects given by the system 



Here, u is a function of the variables (x, t) to be thought of as a density or concen- 
tration, and therefore nonnegative, while p is the pressure, which is related to u via 
a linear operator /C, which we assume to be the inverse of a fractional Laplacian, 
K, = (— A)~ s , < s < 1. The problem is posed for x G M n , n > 1, and t > 0, and we 
give initial conditions 



where uo is a nonnegative and integrable function in R n with compact support or 
fast decay at infinity. Motivation and related ideas for such a problem are given in 
[TT] , where the existence of a weak solution is established for a suitable class of initial 
data that includes all functions Uq with the above properties. Besides, a number of 
the most important properties is proved, like the energy estimates, the bounds in the 
L p spaces and the property of finite propagation that says that compactly supported 
data produce solutions whose support is compact for all positive times. However, the 
uniqueness of weak solutions is a pending open problem but for the case of one space 
dimension, [I]. Comparison theorems, a crucial tool in parabolic equations are only 
available under special circumstances (i.e., for so-called true super- or sub-solutions), 
as shown in [11J . We will also need the results of [9], where we show that solutions 
with integrable data are actually bounded (so-called L 1 -L°° effect) and furthermore 
that bounded solutions are C a continuous with Holder exponent depending on the 
equation and local Holder constant depending on the L°° norm of the solution. 

In this paper we study the asymptotic behaviour of such weak solutions. First, 
we introduce the rescaled flow and then we prove that the stable configurations of 
this modified equation are just the solutions of an obstacle problem with fractional 
Laplacian, as studied by Athanasopoulos, Caffarelli, Salsa and Silvestre, [2] and [TJ. 
This connection, described in Section [31 is quite surprising. In the original vari- 
ables such equilibrium solution translates into a selfsimilar solution that we propose 
to name "fractional Barenblatt solution" after the analogy with the porous medium 
case, though the analytical properties are quite different, see Theorem 13.11 In Sec- 
tions HI [5] we prove convergence of the rescaled flow to equilibrium by using a certain 





(1.2) 
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"Boltzmann entropy" technique; in the original variables this means that the asymp- 
totic behaviour of finite mass, compactly supported solutions of our problem is given 
by the family of fractional Barenblatt solutions constructed from the stationary frac- 
tional obstacle problem. Uniqueness of the asymptotic profiles is proved in Section 
Eland the proof of the main result, Theorem 14.11 is completed in Section [71 A final 
Section El introduces the open question of spectral gap. 

Notation. We will use the notation L s = (— A) s with < s < 1 for the fractional 
powers of the Laplace operator defined by Fourier transform on the Schwartz class of 
smooth functions in M. n and extended in a natural way to functions in the Sobolev 
space H 2s (M. n ). Technical reasons of paper [IT] imply that in one space dimension the 
restriction s < 1/2 will be observed. The inverse operator is denoted by /C s = (— A)~ s 
and can be realized by convolution 

JC S = K s *u, K s {x) = c(n, s)\x\ 2s ~ n . 

1 /2 

see [TT]. JC S is a positive self-adjoint operator. We will write H s = JC S which has 
kernel K s i%. The subscript s will be omitted when s is fixed and known. For functions 
that depend on x and t, convolution is applied for every fixed t with respect to the 
space variables. We then use the abbreviated notation u(t) = u(-,t). 

Related works. Some papers have appeared since 2008 dealing with similar equa- 
tions. As closely related, we mention the recent research note by Biler et al. [3] where 
a more general version of our model is proposed, and existence of solutions, decay 
estimates and selfsimilar solutions are announced. Full details are still to be provided. 
At least in the latter topic, their method is different and interesting. The asymptotic 
behaviour is not tackled. 



2 Existence result and basic estimates 

Definition. We say that u is a weak solution of equation ( II. ip in = M n x (0, T) 
with initial data u G L 1 (R n ) if u e L l (Q T ), K(u) G W%£(Q T ), and uVK{u) G 
L 1 (Qt), and the identity 

(2.1) J J u(<p t - V/C(it) • V0) dxdt + J u (x) (f)(x, 0)dx = 

holds for all continuous test functions <fi in such that V x <p is continuous, and 
has compact support in the space variable and vanishes near t = T. The following 
result is proved in [llj. When T = oo we write Q instead of Qt- 



3 



Proposition 2.1 Let u G L 1 (M n ) n L°°(R n ), w > 0, and swc/i t/iat 

(2.2) u (x) < Ae'^ for some A, a > . 

Then there exists a weak solution u of Equation (II. ip initial data Uq. Besides, 
u G L°°(0, oo : L l (R n )), u G L°°(Q) ; VH(m) G L 2 (Q). For all t > we have 

(2.3) / u(x,t)dx = / u (x)<ia;, 

Jr™ JR™ 

and ||zi(t)||oo < 1 1 Mo 1 1 oo- 27ie solution decays exponentially as \x\ — )■ oo. T7ie /irst 
energy inequality holds in the form 

(2.4) / / |VHw| 2 rfxrft+ / u(ti) log(u(ti)) dx < [ u \og(u )dx, 

JO JR n JM. n ir 

while the second says that for all < t\ < t% < oo 

(2.5) / / u\VK,u\ 2 dxdt + \ I \Hu{t 2 )\ 2 dx <]- f \H{u{t x )\ 2 dx . 

Jti JR n ^ il" ^ </R n 

Properties of the constructed solutions. Here are some of the most useful 

- Conservation of mass 

d f 

(2.6) — / u(x,t)dx = 0. 

dt J 

- LP estimates. We also prove that the L p norm of the solutions, 1 < p < oo, does 
not increase in time. 

- Conservation of sign: uq > implies that u(t) > for all times. 

- Persistence of strict positivity: At every point x where u(xo,to) > for some time 
to we have u(xo, t) > for all later times t > to- 

- Finite propagation. Solutions with compact support: One of the most important 
features of the porous medium equation and other related degenerate parabolic equa- 
tions is the property of finite propagation, whereby compactly supported initial data 
Uo(x) give rise to solutions u(x,t) that have the same property for all positive times, 
i.e., the support of u(-,t) is contained in a ball B R ^(0) for all t > 0. 
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Proposition 2.2 [FT] Assume that u is a bounded solution, < u < L, of equation 
f fTTT]) with JC = (-A)~ s with < s < 1 (0 < s < 1/2 if n = 1), as constructed 
in Theorem 12 . 11 Assume that uq has compact support. Then u(-,t) is compactly 
supported for allt > 0. More precisely, if < s < 1/2 and uq is below the "parabola- 
like" function 

(2.7) U (x) = a(\x\ -6) 2 , 

for some a,b > 0, with support in the ball Bb(0), then there is C(n, s, L,a,T) large 
enough, such that 

(2.8) u{x,t) < a(Ct- (\x\ -b)f 
for x G M and <t <T. 



- A standard comparison result for parabolic equations does not seem to work. This 
is one of the main technical difficulties in the study of this equation. In fact, we find 
special situations where some comparison holds by using so-called true super- and 
subsolutions. 

- Next, we state the L 1 to L°° result that is proved in the forthcoming paper [9]. 

Proposition 2.3 Letu be a weak solution of Problem (II. II) - (11.21) withu G L 1 (R Ar )n 
L°°(M, N ), and Uq decreases exponentially as x — > oo. Then there exists a positive 
constant C such that for every t > 

(2.9) sup \u(x,t)\ < Cr a \\u G f T1 



with a = n/{n + 2 — 2s), 7 = (2 — 2s) /((n + 2 — 2s). The constant C depends only 
on n and s. 

- Regularity. Bounded solutions are C a smooth for some a > 0: The Holder 
exponent depends on n and s and the local Holder constant depends also on the L°° 
norm of the solution. This is proved in [9]. 
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3 Large time behavior. Rescaled flow 



We begin here the study of the large time behavior. As a first step in the analysis 
we will introduce self-similar variables, typical in porous medium theory, leading to a 
rescaled evolution equation. In the analysis of the steady states of that flow we will 
discover the solutions of an elliptic obstacle problem with fractional diffusion. This 
surprising connection is the main novelty of the paper, and the basis for the analysis of 
stabilization to be done in subsequent sections. The precise obstacle result is carefully 
stated in Theorem 13. 1[ 

We take a weak solution u > with integrable and bounded initial data, as con- 
structed in [11]. Inspired by the asymptotics of the standard porous medium equation, 
we define the rescaled flow through the transformation 

(3.1) u{x, t) = (i + ty a v{x (i + ty , t) 

with new time r = log(l + 1). We also put y — x (1 + 1) _/3 as rescaled space variable. 
In order to cancel the factors including t in a explicit way, we get the condition on 
the exponents 

(3.2) a + (2-2s)/3 = 1, 
where we use the homogeneity of K. in the form 

(3.3) ()Cu)(x,t) = (l + t)~ a+2s ^)Cv)(y,T). 

Since we also want conservation of (finite) mass, we must put a = n(3. We get the 
precise value of the exponents: 

(3.4) = 1/(71 + 2-2*), a = n/(n + 2-2s). 

We recall that < s < 1 so that G (l/(n + 2), 1/n). The first value is the standard 
porous medium case. In this way, we arrive at the nonlinear, nonlocal Fokker-Planck 
equation 

(3.5) v r = Vy(v(V y lC(v)+Py)) 

with (3 as given above. Note that this formula implies a transformation for the pressure 

of the form 

(3.6) 

ft — 2s 

p(u)(x, t) = (l+t)- a p(v)(x (1+t)"' 3 , r), with a = a-2s/3 = 1-2/3 = — — < 1. 

t~t I 2. 2iS 
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We are not making much use of this last formula. 

• Equilibrium states for the rescaled flow 

We want to find stationary solutions of the rescaled equation (13.50 . i.e., solutions 
V(y) of the system 

(3.7) V,- (VV y (P + a\y\ 2 )) = 0, P = K{V). 

where a = (3/2, and (3 defined just above. Since we are looking for asymptotic profiles 
of the standard solutions of (II. 10 we also want V > and integrable. The simplest 
possibility is integrating once and getting the radial version 

(3.8) VV y (P + a\y\ 2 )) = 0, P = K(V), V > 0. 

The first equation gives an alternative choice that reminds us of the complementary 
formulation of the obstacle problems, [12], [5]. 

• Obstacle problem. Barenblatt solutions of new type 

Indeed, if we solve the obstacle problem with fractional Laplacian we will obtain a 
unique solution P(y) of the problem: 

r39l P>$, V = (-AYP>0; 

^ ' ' either P = $ or V — 0. 

with < s < 1. In our present application we have to choose as obstacle $ = C—a\y\ 2 , 
where C is any positive constant and a = (3/2. For uniqueness we also need the 
condition P — > as \y\ — > oo. The theory is developed in the papers [21 [7], the 
solution of this obstacle problem is unique and belongs to the space H~ s with pressure 
in H s . The solutions have the following regularity: P E C 1,,s (IR n ) and V E C 1_s (IR n ). 
The figure represents an approximate plot of P{y). 
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Note that for C < the solution is trivial, P = 0, V = 0, hence we choose C > 0. 
We also note that the pressure is defined but for a constant, so that we may take 
without loss of generality C = and take as pressure P = P — C instead of P. But 
then P — > implies that P — > —C as \y\ — > oo, so we get a one parameter family of 
stationary profiles that we denote Vc{y). In any case, if we consider the free boundary 
location, |x| = R(C,a,n), which is the boundary of the contact set of the obstacle 
problem, C = {\y\ < R(C)}, then R tends to as C — » 0, and R tends to oo as 
C — > oo. Let us summarize the results of this section. 



Theorem 3.1 For every C > there exists a unique solution P = Pc{y) £ C 1,s of 
the obstacle problem (13. 9p with $(y) = C — a\y\ 2 and P — > as \y\ — > oo. It is 
radially symmetric as a function of y and we have the scaling law 

(3.10) P c {y)=CP 1 {y/C l l 2 ). 

Moreover, Vc = (—A) s Pc satisfies the scaling relation Vc(y) = C 1 ^ s Vi(y/C 1 ^ 2 ) and 
V\(y) has compact support. Finally, if we put a = (3/2 > 0, then 

(3.11) U c (x,t) = (l + t)- a Vc(x(l + t)-P) 
is a weak solution of Equation (11.11) . 



Note that Uc{x, t) is a weak solution of Equation (11.11) in the sense we have defined. 
It lives for t > — 1 and it takes at the initial time t — —1 a Dirac delta M 5(x) as 
initial data, i.e., it is the source-type or Barenblatt solution for this problem, and it 
has a profile Vc > that has compact support. In the sequel we call these solutions 
the family of fractional Barenblatt solutions. 

Some further observations: an easy calculation gives the realationship 

(3.12) M = Ju c (x,t)dx = Jv c (y,r)dy = c(n,s)C {n+2 - 2s)/2 

(compare with the PME case in [12], page 23). Clearly, V > is supported in the 
contact set of the obstacle problem, C = {\y\ < R(C)}. The radius R is smaller than 
the intersection of the parabola $ with the axis R 1 = (C/o) 1/2 . On the other hand, 
the pressure Pc i\y\) is always positive and decays to zero as \y\ —> oo according to 
fractional potential theory, [151 02] • 
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4 Entropy estimates for the rescaled problem 



Our next problem is now to prove that these profiles are attractors for the rescaled 
flow that we have introduced. In the next sections we prove the main result of this 
paper. 

Theorem 4.1 Letu(x,t) > be a weak solution of Problem fll.ip — fjl.2p with bounded 
and integrable initial data such that uq > has finite entropy in the sense defined in 
formula (I4.4j) . Let v(y,r) be the corresponding rescaled solution. As r — > oo we have 



The constant C is determined by the rule of mass equality: f Rn v(y, r) dy = L n Vc(y) dy. 
In terms of u this translates into 



u(x,t) — Uc(x,t) in L 1 (M n ), t a \u(x, t) — Uc(x, t)\ — > uniformly in x , 
both limits taken as t — > oo. 

As a first step in the proof, we review the estimates of Section [2] in order to adapt 
them to the rescaled problem. There is no problem is reproving mass conservation 
or positivity for the rescaled flow, we leave it to the reader checking both facts. By 
shifting a bit the origin of time, there is no loss of generality in assuming that the 
entropies mentioned below are finite even at t — even if they were not assumed to 
be so. The values of a and (3 are those fixed in (13. 4ft . 

• The first energy estimate becomes (recall the notation "H = /C 1 / 2 ) 



• However, the second energy estimate has an essential change. We need to define 
the entropy of the rescaled flow as 



(4.1) 



v(-,t) -»> V c (y) in 



L 1 (M n ) and also in L°°(R n ). 



(4.2) 



(4.3) 




(4.4) 




The entropy contains two terms. The first is 
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hence positive. The second is the moment £ 2 (v(t)) = M 2 (f(r)) := f \y\ 2 v dy, also 
positive. By differentiation we get 



(4.5) 



-X(v) 



X(v) 



V{lCv + '-\y 



vdy. 



We call X(v) the entropy dissipation of the rescaled solution v. This is formal, the 
rigorous formula says that whenever the initial entropy is finite, then £{v(r)) is uni- 
formly bounded for all r > 0, X(v) is integrable in (0, oo), and 



(4.6) 



S(v(t)) + 



vdy dt < £(vq). 



The bound on £{v{t)) implies a bound for the energy £i{v{r)) and for the moment 
M 2 (v(t)) in L°°(0,oo). We also conclude that 



(4.7) 



dr 



dyv\V(1Cv + (f3/2)\y\ 2 )\ 2 < oo. 



This estimate shows that the integral X(v(t)) converges to zero as r — > in some 
kind of time average, which is a basic fact in the asymptotic analysis. 

Other consequences: (1) For the original variable: the fact that J v(y,T)\y\ 2 dy is 
bounded uniformly in time implies in terms of the original variable u that 



J u{x,t) x 2 dx < Ct 2 K 



(2) The fact that J \H{v)\ 2 dy is bounded uniformly in time implies that 

\U(u)\ 2 dy < Ct 2sfi - a . 



5 Asymptotic behavior. Zero entropy dissipation 

We take a solution u(x,t) under the same assumptions on the initial data, we rescale 
it as in the previous sections to obtain v(y, r), and we put / vo(y) dy = M > 0. Now 
the idea is to let r — > oo in the rescaled flow equation (I3.5p . Since the entropy (14.41) 
is nonnegative and goes down, there is a limit 

E* = lim £{t) > 0. 

t— >oo 
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Moreover, v is bounded in L y uniformly in r, and also J \ VH(v )| 2 dy is uniformly 
bounded, in time, therefore we have that v(-,t) is a compact family and there is a 
subsequence v(-,Tj) with Tj — > oo that converges in L y and almost everywhere to a 
limit v* > 0. 

In this section we prove the following intermediate result: 

Lemma 5.1 the mass ofv* is the same as the mass ofv(-,r), f v^{y) dy = M (hence, 
in particular, is not trivial). The functions v*(y, r) and w(y, r) = /Ct>* + | \y\ 2 are 
continuous and w(-,t) is constant on every connected component of the set where v* 
is not zero, 

Proof. We want to use compactness of the family of time-translates 

Vj(y,r) = v(y,r + Tj), Tj oo. 

We observe that the L°° estimate of Theorem (12.31) with sharp exponents implies 
that v (y, t) is uniformly bounded in time (for all times r > 1), hence the family Vj is 
uniformly bounded. Moreover, bounded families of solutions are uniformly equicon- 
tinuous (Holder continuous), and this applies both to u and v. This was proved as 
part of the local regularity of paper [9]. Therefore, the convergence Vj — > v* takes 
place locally uniformly in Q = R n x (0, oo). 

Next, we prove the conservation of mass. This a consequence of the uniform esti- 
mate of v(-,t) in L 00 (R™) for all large t, the local compactness of the solutions, and 
finally the uniform boundedness of the moment J v(y,r) \y\ 2 dy. Another immediate 
consequence is that the lim inf of the entropy £{v{rj)) is equal or larger that M 2 (f*). 

We also have T-L(v) € L 2 uniformly in r as well as VH{y) in L 2 T . Since JC(v) = 
H{H{v)) and V/C(«) = VH{H{v)) = U{V{Uv)), we derive from the bound for VHv 
in L 2 estimates for K{v). We recall that VH{y) is a "derivative of order 1 — s of v" , 
and since v is bounded, v G L^ T , we conclude that v G L 2 H^ C . By potential theory, 
it is then clear that )C(v) G L\H\+ S . 

All of this can be used in passing to the limit in the term J J viVKv^ipdxdt as 
follows: we have the convergence of Vj in C([0,T] : L 2 (i?/j)) together with the weak 
convergence of pj = )Cvj and Vpj. In this way, we find that v* is a weak solution of 
the rescaled equation (13.51) . 

Finally, we want to prove that this solution is stationary. For that we pass to the 
limit in the entropy dissipation bound that says that for every h > fixed 

2 

Vj(y,r) dydr = 0. 



lim 



vtW^ + ^M 2 ) 



n 



By the strong and weak convergence results that we have, it follows that 

(5.1) T{v.) = J J \V(lCv*(y,T) + ^\y\ 2 )\ 2 v^y,r),dydr = 0. 

This implies that if w(y,r) = K.v* + | \y\ 2 then -y*|Vu>| 2 = a. e. in Q. It follows 
that for almost every time the continuous function w must be a constant in space in 
every connected component of the set where the continuous function v* is not zero, 
w(-,t) = K{j). By continuity of w, this happens indeed for all times, but the space 
constant might still vary continuously in time. 

The determination of the space constant, its independence of time, and finally the 
precise form of v * will allow us the unique identification of the limit as a solution of 
the Barenblatt obstacle problem with the same mass J v*(y) dy = M. This is done in 
the next section, and the proof of the asymptotic theorem is completed in Section [7J 



6 Uniqueness for the asymptotic limit 

We will prove here that the vanishing of the entropy dissipation X characterizes the 
stationary solutions of the rescaled flow. We use the notations u and x instead of v 
and y in this "elliptic section". We also put (3 = 1 for simplicity. Thus, we assume 
that u(x) > is continuous in R n , that p = (—A)~ s u is also continuous, and that 
we have the basic identity 

/|V (? -Q)|^ = 0, 

where Q(x) = —\x\ 2 /2. Then 

Theorem 6.1 There is a Cq > such that p equals P(cq), the solution of the frac- 
tional obstacle problem with obstacle $ = Q + cq. Moreover, P(cq) — > as x — >■ oo. 

We proceed via a series of lemmas. We write P — -P( c o)- 

Lemma 6.2 p attains its maximum only at x = 

Proof. By the maximum principle, p cannot attain its maximum at a point where u = 
0, since then (— A) s p = there, and this is impossible at a maximum of p according 
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to the integral formula that describes the operator (— A) s (unless p is constant, which 
is not the case). Hence, p must attain its maximum at a point where u > 0. Besides, 
at the maximum x we have Vp(x ) = 0, hence VQ(x ) = 0, and this means that 
x = 0. 

Corollary 6.3 (i) We have u(0) > 0. (ii) p = Q + c in a neighborhood £1 of x = 0. 

Lemma 6.4 Let P is the solution of the fractional obstacle problem for data $ = 
Q + c 0; then P — Q is strictly convex. 

Proof. Let e be a unit vector and let h > 0. We have 

\ [P (x + he) + P (x - he)] > \ [2c - \{\x + he\ 2 + \x - he\ 2 ) ] = 
c - ^\x\ 2 - \h 2 = Q(x) + c - |/i 2 . 

Therefore, since P is a supersolution of (— A) s p(x) > 0, also the function 
j(x, h) = ^ [Pq(x + he) + P {x - he)} + hi 2 

is a supersolution, and 7(2;, h) > Q(x) + c . Now, the solution of the obstacle problem 
has the property of being the minimal supersolution lying above the obstacle. Hence, 
l{x, h) > Pq(x). Dividing by h 2 we get D ee P > — 1 for any direction e. Now, D ee P 
is a supersolution outside the coincidence set, so D ee P > —1 there. □ 

Lemma 6.5 On each component of the set {u > 0} we have 

p = Q + c 

for some constant (depending possibly on the component). In fact, it holds on the 
closure of each component. We also have p — > as \x\ — > 00. 

The proof is easy using the previous arguments. 

Lemma 6.6 If P is the corresponding solution of the obstacle problem for data $ = 
Q + Co, then Pq >p and Pq = p in the coincidence set of P . 
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Proof. Let P t = Po + t so that Pt — > t > as \x\ — > oo. Then, P t — Q is strictly convex 
outside the coincidence set (see proof above), thus the gradient never vanishes. The 
obstacle is now $ t = $ + 1 so there is no coincidence between P t and $ = Q(x) + c . 
Let A be the coincidence set of Pq, the closure of a ball B p (0), and let 0o be the 
connected component of {p = Q + cq} near the origin. 

From the previous estimates, for t large Pt does not touch p. Let us now decrease t 
and suppose that there is a first to > such that P to touches p. 

The possibility of first contact P t = p at a point xq such that u(xq) = contradicts 
the maximum principle, unless Pt and p coincide and then to = 0. 

The first contact cannot happen either at a point xq where (— A) s Pq = and u > 
since then p = Q + c\ [c\ > cq) in a neighborhood of xq, and we have already proved 
that VPt cannot equal VQ. 

Finally, we have to consider the case (— A) s Pq > and u > 0, where Pt = Q + cq + 1 
while p = Q + c. This forces again to = 0. This concludes that proof that Pq > p. 

In order to check that Ao C 0o we consider for contradiction the existence of a point 
xo G B p (0) such that xq G OGq. In that case Pq — p has a zero minimum there with 
(-A) s P (x ) > 0, (-A) s p)(x ) = 0, hence (-A) S (P - p)(x ) > which contradicts 
the maximum principle. □ 

Lemma 6.7 We have p>Q + c . 

Proof. Suppose that p — Q — cq has a negative minimum at some point Xq. We can 
translate p up and bring it down until it touches Pq by above. This must happen 
on the set where Uq = 0, since A C 0o- But this contradicts the maximum principle 
since u > at x and Pq and p are not equal. 

Finally, we prove that p is a supersolution and p> P (the least supersolution). □ 



7 End of proof of Theorem 14.1 



We now go back to the end of Section [5] and revert to the notations y and v. Once we 
have uniquely identified the pressure p* at a.e. time as the solution P Q of the obstacle 
problem for a certain constant Cq that might depend on time, the scaling laws (cf. 
( 13.1 2p ) allow to uniquely identify c as the constant such that J (—A s P )(y) dy = M. 
Hence, cq is fixed in time and we conclude by continuity that p*(y, r) = P Co (y). 



14 



The above argument also identifies in a unique way the constant asymptotic limit 
t) = Vc(y)- It follows that the convergence Vj — > V c as Tj — > oo takes place 
independently of the subsequence Tj — > oo, and we conclude that Vj(-,r) — > Vc as 
r — > oo. Since the families Vj are compact in iL^ c (M n ) the convergence is uniform in 
that space. On the other hand, from the conservation of mass and the uniform bound 
for J \y\ 2 v(y,r) dy we conclude that the convergence takes place also in L 1 (M n ). 

We still need to check the uniform convergence of v(-,t) to Vc in the whole space. 
This is a rather simple calculus lemma using the following facts: 

(i) for r > 1 the funtions v(-,t) are uniformly bounded, 

(ii) for every e > there exists R = R(e) such that j^ >R ^ v(y, r) dy < e uniformly 
in t. 

(iii) uniformly bounded functions are C a continuous with a uniform bound. 

This ends the proof for v. The translation of the results in terms of the variable u 
is immediate. 



8 The question of a spectral gap 



When using the entropy-entropy dissipation approach, we are interested in the fol- 
lowing inequality 



8{v) < CI(v) 



for a class C of functions that includes our solutions and for a constant C > which 
is called the spectral gap. In our case it means that for all v G C we must have 



V(/C^ + ||y| 2 ) 



vdy. 



(8.2) f \Hv\ 2 dy + (5 [ \y\ 2 vdy<c[ 

Developing the last integral leads to 

Hv\ 2 dy + [3 [ \y\ 2 vdy<C [ | V/Cw| 2 vdy + C(3 2 I \y\ 2 vdy + 2C/3 I (VJCv.y)vdy. 



We can transform the last integral into 

J (Vfcv.y) vdy = —n J Kvvdy — J (K,v)(Vv.y) dy. 
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For many solutions the last term will be positive, hence we need 

(l + 2Cn/3) J \Uv\ 2 dy<C J \ VKv\ 2 vdy + {Cfi 2 - (3) J \y\ 2 vdy. 

The question is: can we prove this strange interpolation inequality? Since the order 
of differentiation in the first term is — s, for the second 1 — 2s, while the third is zero. 
It is clear that the case s = 1/2 is different and maybe simpler. 
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